Experimental estimates of dephasing time in molecular magnets. 
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Muon spin relaxation measurements in isotropic molecular magnets (MM) with spin value S 
ranging from 7/2 to 27/2 are used to determine the magnitude and origin of dephasing time of 
molecular magnets. It is found that ~ 10 nsec with no S or ligand dependence. This indicates a 
nuclear origin for the stochastic field. Since T4, is a property of the environment, we argue that it is 
a number common to similar types of MM. Therefore, T0 is shorter than the Zener and tunneling 
times of anisotropic MM such as Fes or Mn4 for standard laboratory sweep rates. Our findings call 
for a stochastic Landau-Zener theory in this particular case. 

PACS numbers: 75.50.Xx, 76.75.+i 



Quantum tunneling of the magnetization in anisotropic 
molecular magnets (MM) with high spin value is a fasci- 
nating subject which contrasts clean and accurate experi- 
mental data with sophisticated theoretical models [l| . At 
the heart of these models stands the Landau and Zener 
(LZ) derivation of quantum tunneling between levels, 
which at resonance have a tunnel splitting A, but are 
brought into and out off resonance by a time-dependent 
field. This model can be described by the Hamiltonian 
= ptSz + ASx, where S is the electronic spin opera- 
tor and P is proportional to the external field sweep rate 
dH/dt. The LZ theory predicts the transition amplitude 
Clz that a spin prepared at time < = — 00 in the low en- 
ergy state |-|-) will be in the high energy state at t = 00 
which is again |-|-), namely, Clz = where U 

is the time propagator operator. The calculation of this 
amplitude has a path integral representation as demon- 
strated graphically in the inset of Fig. [H [1, @|. In 
this inset, the solid lines show the instantaneous energies 
E± = ±i-\/A2 -|- p^t'^j and a single path is associated 
with a transition from the lower energy state to the up- 
per energy state, which occurs at a specific time t' . 

However, there is consensus among researchers that 
the tunneling in MM is incoherent due to interactions of 
the spin with a stochastic field B {t) which is produced by 
nuclear moments [i, Si, Bin, and that the dephasing 
time of the quantum states must be taken into account. 
The dephasing time is defined using the corelator of 
the stochastic field 



(B(t)B(0)) = (B2)exp(-i/Tc), 
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probability Plz is given by the absolute value square, of 
the sum of the transition amplitudes, for different paths. 
This yields the famous expression 



Plz = 1 - cxp 
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When the dephasing time is very long, the transition 



of flipping states [ll]. In contrast, if the dephasing time is 
very short, the interference between paths should be de- 
stroyed and transition probability should become a sum 
of instantaneous transition probabilities. 

Therefore, there arc four important time scales in the 
LZ problem: I) the tunneling time It = h/ A which is 
set by the tunnel spliting, II) the Zener time tz = A//3, 
which is the time segment around t = where tunneling 
can occur during a field sweep in the adiabatic case {It ^ 
tz). III) the correlation time Tc, and IV) the dephasing 
time over which different paths interfere coherently. 
Determining these time scales even roughly could help 
select the theory for the analysis of magnetization jump 
experiments. Moreover, theories are available only for 
particular orders of time scales, which might not be the 
realistic ones. 

The theories addressing the stochastic LZ problem 
can be divided into two groups according to the type 
of stochastic field they use: Ising type with coupling 
Bz{t)Sz, or Heisenbergtype with an B(i)-S term. In the 
Ising case Kayanuma found modifications to the LZ 
formula for the order of time scales Tc ^ <C {tztTY^"^ 
and Tc ■cMax[tz, i^iT/T^]. In this case the transition 
probability is give by P = [1 — exp {—Titz/tT)] /2. There- 
fore, when the transition is sudden (t^ 3> tz) then 
P = 7rA^/(2fi,/3) as in Eq. [3]at the same limit. Sinitsyn, 
Prokof 'ev, and Bobrovitski ^ extended this work using 
macroscopic spin bath description of Bz{t) and showed 
that Kayanuma's sudden result is correct if and only if 
tr ^ tz- In the Heisenberg case, Shimshoni and Stern 
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found corrections to the LZ formula in all orders of time 
scales they examined. Here we mention just the interest- 
ing case of Tc ^ T0 and tx ^2 where they found 
that P ~ 1 - (T^/tz) [exp{2tz/T^)PLZ + (tT/T^)^] 0. 
More theoretical work can be found in Ref. hj. The con- 
sensus seems to be that when the field sweep is adiabatic 
tz ^ tr, the stochastic field modifies the LZ formula, 
and that in the sudden limit It 3> tz of the Ising case 
the dephasing time has no impact on the tunneling 
probability. However, as far as we know there is no the- 
ory for the Heisenberg coupling when is the shortest 
time scale in the problem. 

Despite the importance of determination in the LZ 
problem, today there is no experimental estimate of this 
time in the problem of magnetic quantum tunneling. The 
purpose of the present work is to provide such an esti- 
mate. We do so by measuring the dephasing times of 
isotropic molecular magnets (A = 0) with different spin 
value and ligands, and project the result to anisotropic 
MM such as Fcg or Mn4. This allows us to set the order 
of tz, tr, and r^. Our major finding is that is the 
shortest time scale in the problem. Since nuclear dipolar 
coupling to the molecular spins involves all directions, 
we conclude that there is no relevant theory for the LZ 
problem in MM with stochastic field. 

In addition to the contribution to the problem of mag- 
netic quantum tunneling, our experiment has its own 
merit. It is the first examination of magnetic fluctuation 
as a function of the spin value S. As such it provides a 
new look at the interaction between spins and the lattice 
in the quantum (temperature independent) regime. 

We determine the dephasing times of isotropic 
molecules by performing muon spin relaxation measure- 
ments on eight different MM with A ~ and spin value 
ranging from S = 7/2 to S* = 27/2. The major as- 
sumption here is that is a property of the environ- 
ment and not of the molecule (see Eq. [2]). Therefore, if 
we determine for one type of molecule, and if a dif- 
ferent molecule has the same environment, it will have 
the same r^. This assumption received experimental 
support recently in the work of Ardavan et at. They 
showed using ESR that two different molecules, one with 
zero field splitting and the other without it, have the 
same electronic T2 However, it also has advocates. 
Stamp, Tupitsyn, and Morello argue that the molecu- 
lar electronic spin impacts the nuclear spin dynamic and 
therefore the dephasing time should depend on A so that 
T0 oc A Q. Encouraged by the experimental finding we 
continue the presentation using our assumption. 

What allows us to extract the dephasing time is the 
fact that our muons are coupled to the electronic spins 
of isotropic MM that experience only the stochastic and 
external fields. Therefore, the number of parameters 
needed to be determined in our experiment is minimal, 
and there is no need to know a priori the order of time 
scales. The leading terms for such an Hamiltonian are 

?^ = -2/XB[H + B(<)]S + fi^7[H + SA]I (4) 
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FIG. 1: (color online). The muon spin lattice relaxation rate 
l/Ti as a function of temperature and field in the CrNi2 
{S = 7/2) high spin molecular magnet. Inset: the solid lines 
show the instantaneous energy levels as a function of time in 
the Landau-Zener problem. Dashed line is a schematic rep- 
resentation of a path the spin can take when tunneling from 
the low energy state to the high energy state at time t' /tz- 



where I is the muon spin, H is the external field, ^7 = 
851.62 MHz/T is the muon gyromagnetic ratio, /i^ is 
the Bohr magneton, and A is a coupling matrix. We 
ignore the B(t)I term since the field experienced by the 
muon from the molecular spins is greater than this term. 
Due to the fluctuating field B, S will vary in time. The 
simplest assumption that one can make is that the cor- 
relation function ({S(i), S(0)}), where {} stands for an- 
ticommutator, decays exponentially. The decay rate is 
determined by the dynamic properties of B(i) which is 
produced by the environment of the molecules. There- 
fore, we expect 

{S(t),S(0)} = 252exp(-VT^) (5) 

with set by Eq. [51 It is possible that will be H 
dependent but we will show experimentally that this is 
not the case for < 2 kG. 

We investigated CrCu4 (5* = 7/2), CrNia (S = 7/2), 
CrNi2Mn4 {S = 13/2), CrNi2Ni4 {S = 15/2) and 
CrNiMns {S = 20/2). To this, we added data from a 
previous study of CrCug (S = 9/2), CrNig (S = 15/2), 
and CrMng {S = 27/2) by Salman et al. Q. These com- 
pounds, based on polycyanometalated precursors, are 
prepared following a step-by-step synthetic strategy. The 
key idea is to use polydentate amine ligands in order 
to avoid polymerization and get discrete entities with 
well-defined spin and anisotropy [l^, [l^. Most of the 
compounds are fully described in the literature [13, [3 ■ 
They may be divided into two groups: i) isotropic high 
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spin molecules (CrCug, CrNig, and CrMng ), ii) nearly 
isotropic molecules with no detectable energy gap or 
smah one ~1 K (CrCu4, CrNia, CrNiMns, CrNi2Mn4, 
CrNi2Ni4). 

In our fiSK-Ti experiments we measure the polariza- 
tion P{t, H) of a muon spin implanted in the sample, as a 
function of time t and magnetic field H, when the field is 
applied in the direction of the initial muon polarization. 
These experiments were performed at both ISIS and PSI, 
exploiting the long time window in the first facility for 
the slow relaxation of the low S molecules, and the high 
time resolution in the second facility for the fast relax- 
ation of the high S molecules. Typical raw /xSR data are 
presented in Ref. [13] ■ The data for all samples are fitted 
to P{H, t) = exp(— -^/t/Ti) -I- Bg, where Bg is a field and 
temperature independent background. This root expo- 
nential behavior is a consequence of the many different 
muon sites in the sample. 

In Fig. [T] we depict the temperature dependence of 
l/Ti. As the temperature is lowered, the relaxation in- 
creases due to slowing down of the spin fluctuation as a 
result of the interactions between spins in the molecules. 
However, once the MM is formed, the spin dynamics is 
nearly temperature independent down to the milikelvin 
regime. All molecules show the same behavior. More raw 
Ti data are presented in Ref. 14]. 

In Fig. [2] we depict Ti as a function of for all the 
molecules measured to date. There is a large variation 
in the scale of Ti between the different molecules. A 
linear dependence of the form Ti = m + nH^ is found in 
all cases, as demonstrated by the fitted solid line. This 
is in agreement with Ref. [1^. The difference between 
molecules is in the slope n and crossing of the line m. 
The dephasing time could be extracted from the standard 
theory of Ti relaxation where 
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Although this expression is a result of perturbation ex- 
pansion where H provides the quantization axis, it was 
demonstrated by numerical methods to be a good ap- 
proximation even for 7? — > ^19.] . Here we assumed for 
simplicity that A is diagonal and isotropic, but this as- 
sumption has no significance for our conclusions, is 
obtained from 



T4> 



1/2 



(7) 



for each molecule at the lowest temperature. 

The main experimental finding of this work is pre- 
sented in Fig. [3] where t^^ is plotted as a function of 
S for all the molecules. This plot shows that within ex- 
perimental errors is weakly dependent on the type of 
molecule used, despite the large variations in Ti. In par- 
ticular, is weakly dependent on S or the ligand. To 
emphasize this conclusion we fit the data to three differ- 
ent power laws: t7^ oc constant, S, and S"^. The quality 
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FIG. 2: (color online). Muon Ti at 100 mK versus field 
squared for all molecules including three from Ref. p^ . The 
solid lines are linear fits. 



of the fit expressed as the value of the reduced is 
shown on the graph. The = const gives an order 
of magnitude better fit than the other power laws. It is 
also interesting to compare our finding of ^ 10 nsec 
to other experiments. In the deuterated molecules CryMi 
and CryMn the ESR T2 (interpreted here as r^) is 3 //sec 
p^ . Had the samples were not deuterated, would have 
been 80 nsec due to the gyromagnetic ratio between pro- 
tons and deuterium. In the V15 molecule ~ 2 nsec 

M- 

It is highly significant that is nearly spin- and 
ligand-independent. Since is determined by the en- 
vironment in which the molecules are embedded, its 5- 
independence means that coupling to other molecules or 
to phonons is not responsible for t^. In both these mech- 
anisms the field B experienced by a given molecule de- 
pends on 5", and according to Eq. [2]we would expect an 
5-dependent r^. We therefore conclude that at T ^ 
the stochastic field B(i) responsible for the MM spin mo- 
tion emanates from nuclear moments, most likely pro- 
tons. Since there are many protons in the ligands, the 
variations between ligands do not have a big impact on 
T0. According to Eq. [21 on the order of 10 nsec could 
be generated by a field B ~ 1 to 0.01 G, which for 5 = 10 
is equivalent to 200 to 0.2 MHz, fluctuating at a rate of 
1/tc ~ 4 to 4 X 10~^ /xsec"-^, respectively. These values 
are typical for nuclei. In Feg Morello et al. found nuclear 
I/T2 on the order of 10"^ /tsec"! [2l|. 

As we argued before, the dephasing time should be 
typical of high spin magnetic molecules made of transi- 
tion metal ions embedded in a sea of protons. Indeed, 
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FIG. 3: (color online). Dephasing rate t^^ extracted from 
the muon relaxation data as a function of spin value S of the 
various magnetic molecules. The solid lines are fits to power 
laws as indicated in the figure, x represent the quality of the 
fit. 



the eight isotropic molecules reported here are different 
but have similar r^. We have no experimental reason to 
believe that will be substantially different in Fes or 



Mn4 where A was measured. In both cases A ~ 10"'' K 
for the —S to S transitions The tunneling time 

tx — h/A ~ 7.6 X 10~^ sec. This tunneling time is 
longer than the dephasing time ~ 10~® sec in our, 
and other [l^, [23| molecules. Moreover if, for example, 
(3 = 0.001 K/sec, then = A/P - 1 x 10^-* sec (for the 
same transition). This implies the order of time scales 
tz tx ^ T^T St regime which corresponds to a strong 
dephasing. For /3 = 0.1 K/sec we have tx ^ tz S> r^. 
As we mentioned before, the impact of the Heisenberg 
type stochastic fluctuations in this order of time scales 
on transition probabilities is not known theoretically. 

To summarize, we have measured spin correlations in 
isotropic molecular magnets on a wide range of S val- 
ues. We found that the correlation time is nearly S- 
and ligand-independent and on the order of 10 nsec. We 
use this time as an estimate of dephasing times in non 
isotropic molecules such as Fcg and Mn4 where tunneling 
occurs. Our findings call for a theoretical development of 
the LZ problem with stochastic field fluctuations coupled 
to all components of the spin S operator, where is the 
shortest time scale in the problem. 
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